Single spin asymmetries in hadron-hadron collisions 
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We study weighted azimuthal single spin asymmetries in hadron-hadron scattering using the 
diagrammatic approach at leading order and assuming factorization. The effects of the intrinsic 
transverse momenta of the partons are taken into account. We show that the way in which T-odd 
' functions, such as the Sivers function, appear in these processes does not merely involve a sign 

flip when compared with semi-inclusive deep inelastic scattering, such as in the case of the Drell- 
Yan process. Expressions for the weighted scattering cross sections in terms of distribution and 
fragmentation functions folded with hard cross sections are obtained by introducing modified hard 
cross sections, referred to as gluonic pole cross sections. 
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> \ I. INTRODUCTION 

00 . 

^v^j ■ Accessing the effects arising from the transverse momentum of quarks in hadrons requires hard processes involving, 
' at least, two hadrons (or hadronic jets) and a hard scale to separate them. This is most cleanly achieved in electrowcak 
processes in which the gauge boson provides the hard scale separating the two hadronic regions. The transition 
between the hadronic regions and the hard subprocess is described by soft quark and gluon correlation functions, 
implying approximate collinearity between the quarks, gluons and hadrons involved. Without effects of quark intrinsic 
■ transverse momentum these are bilocal, lightlike separated, matrix elements where collinear gluons provide the gauge- 
fi i| link. Transverse momentum dependent correlation functions involve bilocal matrix elements off the lightcone fil. Here 
the issue of color gauge-invariance is slightly more complex, involving gauge fields at lightcone infinity 0, IH, M IE El • 
' The gauge-link structure may depend on the hard subprocess and leads to observable consequences. 

Absorbing the soft physics in the correlation functions requires coupling of, essentially, collinear quarks (and gluons) 
to the hadronic region. These partons themselves are approximately on mass-shell. In the absence of a hard scale 
from an electroweak boson, as in strong interaction processes, the simplest hard subprocess (large momentum transfer) 
involving on-shell quarks and gluons is a two-to-two process. 

In this paper we discuss hard hadron-hadron scattering processes using the diagrammatic approach rather than 
the commonly used helicity approach 01 B S liil liil LL=f • This has the advantage that we can directly connect to 
the matrix elements of quark and gluon fields, without having to go through the step of rewriting them into parton 
distributions with specific hclicities. It allows us to include the effects of collinear gluons, determining the gauge-link 
structure and to compare this for semi- inclusive deep inelastic scattering (SIDIS) and the Drell-Yan process (DY). 
We note that in this paper we assume the validity of factorization. 

We will consider the possibilities to measure transverse moments which are obtained from transverse momentum 
dependent (TMD) distribution and fragmentation functions upon integration over intrinsic transverse momentum 
(k T ) including a fc T -weighting. In the transverse moments the effects of the gauge-link structure remain visible. For 
this one needs to classify the distribution and fragmentation functions as T-even or T-odd. In single-spin asymmetries 
at least one (in general an odd number of) T-odd function appears, while in unpolarized processes or double-spin 
asymmetries an even number of T-odd functions must appear. The importance of considering transverse momentum 
dependence comes from the fact that for spin and spin i hadrons the simple transverse momentum integrated 
distribution and fragmentation functions, relevant at leading order, are all T-even. 

The specific hadronic process that we will consider is the 2-particle inclusive process H1+H2 — > hi+h,2+X, which 
in order to separate the hadronic regions requires minimally a two-to-two hard subprocess. Also included are inclusive 
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FIG. 1: The leading order contribution to the cross section of H1+H2 — > hx+hz+X. 



hadron-jet and jet-jet production in hadron-hadron scattering. The 1-particle inclusive process +p — » 7T+X involving 
a transversely polarized proton is known to show a large single-spin asymmetry [l^, Hj, HU UM LLdl ■ Some of the 
mechanisms Uj& |2(J, Wi UlM to explain these asymmetries involve T-odd functions, such as the Sivers distribution 
function or the Collins fragmentation function [23. These functions are expected to appear in a cleaner way in 
2-particle inclusive processes |lflj . Here we only consider single-hadron fragmentation functions, in which case the 
2-particle inclusive production requires hi and /12 to belong to different, in the perpendicular plane approximately 
opposite, jets. 

In this paper we limit ourselves to the (anti)quark contributions with as main goal to show the relevant gauge-link 

structure for the T-odd Sivers distribution functions f^p and the Collins fragmentation functions H^ 1 ^ entering these 
processes. This is important for the study of universality of these functions. The paper is structured as follows. In 
section HT1 we consider the kinematics particular to 2— >2 particle scattering. In section ITTT1 we discuss our approach and 
several (weighted) scattering cross section are written down for hadronic pion production and hadronic jet production 
in section HVI Details about the gauge- links and their consequences for distribution and fragmentation functions are 
dealt with in the appendices. 



II. KINEMATICS 



?cm 



The hard scale in the process Hi (Pi) +H2 (P2) —> hi{Ki)+h,2{Ki)-\-X is set by the center-of-mass energy *Js = E c 
The leading order contribution to the scattering cross section is shown in Fig. ^ In a hard scattering process it is 
important to get as much information about the partonic momenta as possible, in our case including, in particular, 
their transverse momenta. The partonic momenta, for which p\-P\ ~ p\ ~ Pf = M\ are of hadronic scale, are 
expanded as follows 

pi = X1P1+ (Ji rii + p lT , (la) 
P2 = X2P2+ <?2 n 2 + P2t , (lb) 
ki = z-f 1 A'i + a[ n' x + k\ T , (lc) 
k 2 = z^ 1 K 2 + <t' 2 n' 2 + k 2T , (Id) 

where the rij (nj) are lightlike vectors chosen such that P\-n\ cx 0(s 1 ^ 2 ) and similarly for the other partonic momenta. 
The fractions X{ = prfiij 'Pi-rii and z^ 1 = ki-n^j ' Ki-n^ are lightcone momentum fractions. The quantities multiplying 
the vectors Hi are of order s^ 1 / 2 and are the lightcone components conjugate to p,-n,-. They are given by 

ai - p. ni ' (2) 

with similar expressions for the a[. If any of the 'parton' momenta is actually an external momentum (for leptons or 
when describing jets) the momentum fractions become unity and the transverse momenta and <ii vanish. 
Integration over parton momenta is written as 

d 4 px = dxi d 2 p lT dipx-Px) , (3) 

with d(p\-Pi) = (Pi-ni) dai and similar expressions for d 4 p2, d 4 ki and c? 4 fc 2 . The integrations over the parton mo- 
mentum components (p%-P%) and [kvKi] will be included in the definitions of the TMD distribution and fragmentation 
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FIG. 2: Plane perpendicular to the incoming hadronic momenta. 



functions. Note that the subscripts T have a different meaning in each of the above decompositions, i.e., p± T is trans- 
verse to Pi and rii, while P2 T is transverse to P2 and n 2 , etc. Momentum conservation relates the partonic momenta: 



Pi + P2 - h - k 2 = . (4) 

This relation, however, does not imply that the sum of the intrinsic transverse momenta, q T = pi T -\-p 2T —ki T —k 2T 
« z-f 1 Ki+z^ 1 K2—X1P1—X2P2 vanishes. 

We use the incoming momenta Pi and P2 to define perpendicular momenta Ki± orthogonal to the incoming 
hadronic momenta, ifjX'Pi = Kix-P 2 = 0. For the perpendicular momenta it is convenient to scale the variables 
using Xi± = 2\Ki±\/-^s. For the outgoing hadrons we use the pseudo-rapidities r\i defined by r\i = — lntan(i#i) where 
the 8i are the polar angles of these hadrons in the center-of-mass frame. All invariants involving the external momenta 
can be expressed in terms of these variables: 

P 1 -K x = i^ue-" 1 , P 2 -K 1 = \sx 1± e +r * , (5a) 
Pi-K 2 = \sx21_ e- 1 ' 2 , P 2 ■ K 2 = is x 2 i_ e+" 2 , (5b) 

and PyP 2 = \s. These identities are valid up to subleading order in ^s. The remaining invariant K\-K 2 is not 
independent of the others. To leading order, one has 

Ki ■ K 2 = \sx1_1_X2_L cosh 2 [1(771-772)] ■ (5c) 

To subleading order the outgoing hadronic momenta can now be written as 

{K V P 2 ) Pi + {KvPi) P2 , „ 1 / +„ 1D , - mu \ _l v (a \ 

Ki = v h K 1± = ±x 1± (e +r,1 Pi + e m P 2 ) + K 1± , (6a) 



PrP 
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^ = {K2-P2)P 1 + (K 2 -P 1 )P2 + K ^ = , ( £+V2pi + er , 2p2 )+R2 ^ (6b) 

Pl-"2 

The two perpendicular vectors Ki± and K 2 ± are approximately back-to-back (see Fig.EJ). Sometimes the Feynman 
variables XiF = Kf™ / ^™ 1 ( max ) = Xi± sinh r/i are used in the literature. Another useful variable in writing down 
cross sections is the quantity y, which is defined via the Mandelstam variables of the partonic subprocess, and can be 
related to the pseudo-rapidities, 



U= J_ = (PvK 1 )(P2-K 2 ) _ 1 

s V {PvPi){K v K 2 ) e (w-»») + i- {> 

Dividing Ki± and K 2 ±_ by the momentum fractions one immediately sees from the decompositions of the partonic 
momenta that the vector 

r± = 1 , (8) 

z\ z 2 
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only involves transverse momenta of partons. It is just the small projection of the transverse momentum in the 
perpendicular plane, rj_ rs q T j_. The vectors Kix themselves are not 'small' vectors. They are spacelike vectors with 
invariant length of 0(\fs). In the analysis of the kinematics in the transverse plane the momentum fractions are not 
direct observables. In particular, experimentally it is more convenient to work with the directions of the vectors Ki± 
and the corresponding orthogonal directions (see Fig. [2J 

" ^ e1 N = -- e - T —— = e^e lXv , (9) 



L± \Kix\ ' LN s \K n 

and similarly for K~2A_. As illustrated in Fig. [3 the direction e2i_ is, up to a (small)angle 5(f) = </>2— 0i— 7r oc 0(\p T \/ \fs), 
opposite to eu_. 

The momentum conservation relation in Eq. Q is enforced by a delta function in the scattering cross section. The 
delta function can be decomposed using the basis constructed in the previous paragraph. For R = P1+P2 — k\— k 2 this 
decomposition reads 

S 4 (R) = ±s5(R-P 1 )5(R-P 2 )5 2 (R 1 _) = \s 5{R-P X ) 5(R-P 2 ) 5 2 {q T± -r ± ) . (10) 

The arguments of the first two delta functions involve large momenta and can be used to relate the momentum 
fractions x\ and X2 to kinematical observables. For the latter two delta functions the treatment depends on the 
situation. Using the orthogonal vectors eu_ and e\ N we get, up to 0(1/ \fs), 

R-eix = eix-q T + — , (11a) 

\ z\ z 2 / 2 

-R-euv = e lN -q T - — sm(6<f)) . (lib) 
z 2 2 

In the case of two-hadron or hadron-jet production (z2 = 1), the first delta function implies that at leading order 
x\_l/z\ ~ X2i_/z2 = x±, which is interpreted as the scaled parton perpendicular momentum, x±_ — 2 \ k2±\/\/s. Using 
the variable x± as an integration variable we can write 

S 4 (Pi+P2-ki-k 2 ) = \ / dx± S(xx - ^x^+e^)) S(x2 ~ ^x ± (e-^+e-^)) 

S 2 XiJ_ x 2 ±_ J z z q 2 ) 

xf -i Xl - \ x( -i X -L \ x f e iN-q T xj_ . \ 1 ' 
x S(z 1 )8[z 2 )S[ — -= sin {8<t> , 

which shows that in one- or two-particle inclusive processes we are always left with a convolution of distribution 
and fragmentation functions over one momentum fraction or, equivalently, over the parton perpendicular momentum 
variable x±. The last delta function shows explicitly that sin(<5</>) oc 1/t/s and that it can be used to construct cross 
sections weighted with one component of the intrinsic transverse momentum, i.e. e\ N -q T . 

In the case that K\ = k\ and K~2 — fc2 (he. Zi = 1, ki T — 0), such as in production of a lepton pair in Drell- 
Yan scattering or the (idealized) production of two jets, the delta function 5(R-e\x.) also relates intrinsic transverse 
momenta to observed momenta. Therefore, one can construct azimuthal asymmetries involving two components of 
q T . In fact, the product of delta functions 5 2 (R±) = 5(R-e\±) 6(R-ei N ) can be used to weigh with the transverse 
momenta Pi t +P2t, as they relate q T = Pi t +P2t to q = fci+fc2 in the orthogonal plane. With the natural choice of the 
n-vectors in the case that only two hadrons are involved such that P± T = P2 T = 0, one obtains the familiar relation 
Pit+P2t = q T = q-xiPi-x 2 P2, leading to 

S 4 (Pi+P2-ki-k 2 ) = - a i(xi - p^-) S(x 2 - P^-) 5 2 (p lT + P 2 T ~q T (q,PuP2)) ■ (13) 

III. CROSS SECTIONS 

The scattering cross section for p±p 2 — > h\\i2X (see Fig.^l at tree- level is written as 

1 , , „„ (PKx d 3 K 2 , N 
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where the matrix element is expressed in terms of hard amplitudes and correlation functions (see appendix El • It is 
given by 

\A4\ 2 — I dxiSpiT dx2<l 2 p2 T dz^cPkiT dz2 1 d 2 k2 T (27r) 4 <5 4 (pi+p 2 — k\ — fe) 

J ' ' (15) 

x Tr{ ^(xi, pi T )^(x 2 , P2t)^(zi, ki T )^(z 2 , k2 T )H(pi,p2,ki, k 2 )H* (pi,p 2 ,h, k 2 ) } ■ 

The trace involves the appropriate contraction of Dirac indices in soft and hard scattering parts. A summation over 
color and quark flavors is understood. The phase-space elements are given by 

d 3 Ki Xi±s dcj>i 

(2*)*2E Kt =8(2^ ^ Vi aF" (16) 

Combining the phase space integration and the delta functions coming from partonic energy-momentum conservation, 
one has for back-to-back hadron-hadron production 

*r , ( ^ 2 r, 1 „^ 2 «„„^ 2 fc,„^ 2 fc„„ fit 6lN ' qT - fi sinlX^ 

(17) 



da[hih 2 ] = — ^— dx\x dx 2 j_ drjx dr\2 ~P -P f dx± [ d 2 p lT d 2 p 2T d 2 k XT d 2 k 2T s( eiN ^ T - £i S m(S<j)) 
6Zs Zir Zir J J \ y/s I 

x Tr {®(xi,p 1t )$(x 2 ,P2t)A(zi, k lT )A(z 2l k 2T ) H(p 1 ,p 2 , ki,k 2 )H*(pi,p 2 , fci, k 2 )} 



In this expression the momentum fractions are fixed by the arguments of the delta functions in Eq. (|12() . i.e. one has 

Xx (x j_, 771,772)1 x 2 (x±. Tlx, r] 2 ), zi(x±,Xi±), and z 2 (xj_, ^2J_)- Since Xi± < x± and x 2 ±_ < the integration over x± 
is bounded from below. 

In the hadron-jet inclusive process one has A(22,^2t) = 5(z 2 — 1) S 2 (k2 T ) $2 = x± 5(x2±— x±) S 2 (k2 T ) $2, which 
implies that z 2 = 1 and x± — X2± ■ The cross section becomes 

da[hij 2 ] = dxu_ dx 2 ±_ drji drj 2 — — [ d 2 p lT d 2 p 2T d 2 k lT S ( eiN ^ T _ s in(54>) \ 

61 8 Zir Zir J \ y s Z / (18) 

x Tr {$(a;i,pi T )$(x2,P2r)A(zi, k lT ) H(p 1 ,p 2 , k 1 ,k 2 )H*(pi,p 2 , ki, k 2 )} . 

As stated in the previous section, in back-to-back jet production both delta functions in the perpendicular plane relate 
observed kinematical variables to intrinsic transverse momenta. Therefore, in jet-jet production we use the expression 
in Eq. rather than Eq. H12JI for the momentum conserving delta function, leading to 



(19) 



(20) 



da[jij 2 ] = Xl± X2± dxxx dx 2 i_ dr)i dr) 2 ^p- tt^ / ^ 2 p\ T d 2 p 2T S 2 {piT+P2 T -q T ) 
o4 Zir Zir J 

x Tr {G>(x 1 ,p lT )®(x 2 ,P2 T ) H(pi,p 2l fci, k 2 )H*(p 1 ,p 2 , fci, k 2 )} , 

where Z\ = Z2 = 1, X\x = %2± = x± and q T = q—XiPi—X2P 2 - 

In averaged and weighted cross sections we will encounter contractions of hard and soft pieces like 

T,(x 1 ,X2,z 1 ,z 2 , y) = J d 2 p lT d 2 p 2T d 2 k lT d 2 k 2T Tr {$(xi,pi T )$(x2,p2 T )A(zi, ki T )A(z 2 , k 2T ) H H* } 
= Tr {<P(x 1 Mx2)A(z 1 )A(z 2 ) H H* } , 

and 

Eg(xi,X2,zi,Z2,y) = J d 2 p lT d 2 p2 T d 2 k lT d 2 k,2T q" Tr {®(xi,Pit)®{x2,P2t)A{zi, ki T )A(z 2 , k 2T ) H H* } 

= Tr { [ $g(x 1 )$(a; 2 )A(z 1 )A(z 2 ) + <S>{x l )9%{x 2 )A{z 1 )A{z 2 ) 

- $(n)*(a: 2 )Ag(«i)A(z 2 ) - $(n)$(a: 2 )A(2i)Ag(zi)] HE* } . (21) 

These expressions for hadron-hadron scattering (and similar ones for hadron-jet and jet-jet scattering) are schematic 
in the sense that the tracing depends on the particular term in the sum of squared amplitudes, including both direct 
and interference diagrams when the hard amplitude contains more than one contribution. 



(a) (b) (c) (d) 

FIG. 3: Hard scattering amplitudes for quark-quark scattering: (a) t-diagram, (b) w-diagram; 
quark-antiquark scattering: (c) s-diagram, (d) t-diagram. 

In the case of hadron-hadron or hadron-jet cross sections one finds averaged cross sections like 
(da[hiha])= 

J dfo 

dxi± dx 2± drji dr\ 2 #1 f dx± 
= 32^ "2^7 ^7 n^,zuz 2 ,y) , (22a) 

(±sm(8cf>)da[hih 2 ]} = J dfo | s in(^)^^i 

dp_ j dx± eiN ,-£ g ( Xl ^ X2 ^ Zl ^ z ^ y ^ _ (22b) 
Zir J x- ± 



dxi± dx 2 ± dr)i dr] 2 



32tts 3 / 2 



We would like to note that in Eq. (|22b|) . one is weighing with a dimensionless quantity which leads to a suppression 
with 1/y/s. For jet-jet cross sections one has, in principle, the possibility to access both perpendicular directions of 
Eg, assuming that q = k\+k 2 is known accurately. One could, then, weigh with q", in analogy to the Drell-Yan 
process |24j . In that case one weighs with dimensionful quantities, even if these are small momenta, and one does not 
get additional suppression involving the hard scale. 

These equations will be the starting point in the calculation of cross sections. One needs to calculate the quantities 
in Eqs. i|2U[l and 121|) . These expressions involve hard scattering amplitudes and soft correlators <E> and A, which 
are obtained as Fourier transforms of matrix elements of nonlocal combinations of quark and gluon fields. They 
are parametrized in terms of distribution and fragmentation functions as presented in appendices [B] and 10 I n 
order to render the correlators color gauge invariant a gauge-link connecting the fields is needed. In the diagrammatic 
calculation gauge- links are explicitly found by taking into account, for each of the hadrons, the interactions of collinear 
gluons (polarizations along hadron momentum) between the soft and hard parts. Th ese g ive the well-known straight- 
line gauge- links along the lightcone for transverse momentum integrated correlators |25j, but they lead to nontrivial 
gauge-link paths for the TMD correlators 0, 01, H(| . The integration paths in the gauge- links U are process dependent, 
depending in particular on the hard partonic subprocesses. We indicate this dependence by a superscript <S>^ u \x,p T ). 

The transverse momentum integrated correlator is a lightcone correlator with a unique gauge-link, in which the 
path dependence disappears: 

$M(ar) = / d 2 p T $ [u] {x,p T ) = *(a?) . (23) 



For the transverse moments of the correlators obtained after p T -weighting, of which we only consider the simplest 
one, one finds two types of lightcone correlators, a quark-quark matrix element <!>g and a gluonic pole matrix element 
<I>g, where the latter is multiplied by a factor that depends on the gauge-link 

= / d 2 p T V% <S> [U] {x,Pt) = + CP ^{x,x) . (24) 



The gluonic pole matrix element, which contains the T-odd distribution functions, was suggested in a slightly different 
context by Qiu and Sterman H3,|23 as the origin of single spin asymmetries. In processes like SIDIS with underlying 
hard process £+q — > £+q and the DY process with underlying hard process q+q — * £+£, different gauge-link paths 

W'+l andWH appear. In these processes the corresponding factors in Eq. (|24(l are simply Cq ' = ±1. 

As was shown in Ref. |26|. more complex integration paths enter in the gauge- links when other subprocesses are 
involved, such as the two-to-two (anti)quark subprocesses in this paper. Moreover, in general several diagrams enter 
in the calculation. For instance, for quark-antiquark scattering both t- and s-channel amplitudes (see Fig. |3J) can 
contribute, H — H^+H^, whereas for quark-quark scattering we get t- and u-channel amplitudes, H qq — H l qq +H qq . 
In addition, one has to consider the diagrams that produce the gauge-links needed to render the correlation functions 
color gauge invariant. In these cases the gauge- links, in general, also differ for the various terms appearing in the 
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squared amplitude HH^ for a given partonic subprocess. For instance, in the scattering of two identical quarks the 
squared amplitude contains the terms H qq H qq , H qq H q ^, H qq H q £ and H qq H qq . Details are explained in appendix 1X1 
The results of the diagrammatic calculation for transverse momentum integrated cross sections involving soft and 
hard parts can, in leading order, be recast in the form of a folding of the quark distribution and fragmentation functions 
appearing in the transverse momentum integrated correlators <f>(x) and A(z) with hard partonic cross sections. For 
SIDIS one has a folding with the cross section for the hard process £+q — > £+q and in the DY process a folding with 
the cross section for q+q — > £+£. In hadron-hadron scattering one has several hard processes. An example is qq 
scattering with a cross section, in the case of identical quark flavors, of the form 

J A ,«[D] 

aa qq^qq _ OOgq-^qq ,~r:\ 

dt ~ ^ dt ' [ 1 

where the summation is over the different direct and interference contributions involving the t- and u-channel ampli- 
tudes. 

A folding of distribution and fragmentation functions is also possible for weighted cross sections. The cross sections 
involving the link-independent parts of the transverse moments (i.e. $a(x) and Ag(x)) also lead to a folding with the 
normal partonic cross sections, just as for the integrated correlators $(x) and A(z). However, for the contractions 
with the gluonic pole matrix elements it&g and ttAg the gauge-link dependence in the decomposition in Eq. 1241) has 
important ramifications. Expressing the asymmetries as a folding of universal, one argument functions and a hard 
part requires a modification of the hard partonic cross section by including the gauge-link dependent factors Cq^ in 
the various terms in this cross section. This is a convenient way of doing since the values of the factors Cq depend on 
these terms. For instance, in the example of unpolarized qq scattering for identical flavors, the functions appearing 
in the parametrization of the gluonic pole matrix elements are folded with the gluonic pole cross section 

da Tqq _ qq _ C [U(D)] d&qqLqq 

dt ^ G di 

The notation gq emphasizes which quark field (in this case the first one) is accompanied by a zero momentum gluon 
field in the correlator h^>q. The parametrization of this correlator involves one-argument distribution functions, 
which will appear folded with the gluonic pole cross sections. At tree level often only one diagram contributes to 
the partonic cross section. In that case the gluonic pole cross section is simply proportional to the normal partonic 
cross section. For instance, the sign difference between SIDIS and DY for the Sivers distribution function, a uniquely 

defined function appearing in the parametrization of 7t$g(x, x), comes from the factors Cq ' = ±1 discussed above. 
Instead of the folding with partonic cross sections, the Sivers function is folded with the gluonic pole cross sections 

d&eaa^>f.a d& >„_»<>„ 



dt dt 



(27a) 



dt dt { ' 

Although the gluonic pole cross sections should not be interpreted as true partonic cross sections, their concept is 
convenient in order to get a simple folding expression involving the one-argument functions appearing in the gluonic 
pole matrix elements 7t$g and ttAq. Moreover, they are easily obtained from the terms in the hard partonic cross 
section without the inclusion of collincar gluon interactions between the hard and soft parts. 

In the following sections the formalism described above is applied to single spin asymmetries in inclusive two-hadron 
production, hadron-jet and jet-jet production in p^p scattering, for which the gluonic pole cross sections for polarized 
(anti)quark scattering are also needed. 



IV. SINGLE-SPIN ASYMMETRIES IN INCLUSIVE HADRON-HADRON SCATTERING 

As a reference we first consider the cross section for 2-particle inclusive hadron-hadron scattering. The explicit 
expression for the cross section in terms of the distribution and fragmentation functions can be obtained by inserting 
the parametrizations of the correlators, Eqs. I|B7(1 and (|C8|) . into Eq. (|22a|) and performing the required traces. In this 
paper we restrict ourselves to the quark and antiquark scattering contributions. Since the short-distance scattering 
subprocesses remain unobserved, all partonic subprocesses that could contribute have to be taken into account. This 
includes, for a realistic description, besides the (anti)quark contributions qq^qq, qq^qq and 55— >qq, also contributions 
involving gluons, qg^qg, qg^qg, gg^gg, qq^gg and gg^qq including their polarizations |9l flOL HiL H2l 1281 123 . |30 | . 
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However, the (anti)quark contributions suffice to illustrate how the inclusion of gauge-links leads to altered strengths 
of specific distribution or fragmentation functions. Contributions involving gluons can simply be added incoherently 
to the results presented here. For the (anti)quark contributions to the averaged cross section one obtains 



dx± 
x± 



2?r al 

x ; 

9s 



( da[hih 2 ] ) = dx±±_ dx 2 ±_ dt]i dr] 2 — — / ■ 

[ ((l-y) 2 +y 2 )J2ff(x 1 )ff(x 2 )Dl'(z 1 )Dl'(z 2 ) (28a) 
q,q' 

+ (1 ~ % + 1 E n^i)fA^)D\( Zl )Di{z 2 ) (28b) 

+ (1 ~ !,r + 1 E n^i)n{^)D\( Zl )Di{ Z2 ) ( 28c ) 



V 

q,q 

\2 



+ l^hr~ S fK^)fMDl{ Zl )Dl{z 2 ) (28d) 
JirTTEA'^OA'feW^iWW (28e) 
+ ( quark PDFs/FFs antiquark PDFs/FFs ) j + ( K x ^> K 2 ) 

where the summation is over all quark flavors, including the case that q = q' . In this expression y is given by Eq. (J7J 
and s is s = x\ s cosh 2 [ 5(771—772)] = 2^ s/4y(l— y). This result can be recast into a folding of the distribution 
and fragmentation functions appearing in $(x) and A(z) and the elementary (anti)quark cross sections given in 
appendix IEI That is, expression (|28(l can be rewritten to 

(Mhih 2 ])=dx lx dx 2X dr h dr l2 ^ f ^ V /f (xi)/f (x 2 ) | d <^^ Df{z 1 )Df{z 2 ) , (29) 

27T J x± f-^ 2 dt 



where the summation is over all quark and antiquark flavors. In the expressions above the momentum fractions are 
fixed by X\i 2 — ^x±(^e ±Vl +e ±V2 ) and = Xi±/x±. 



A. Hadron-hadron production in p T p scattering: — ► 7T+7T+X 

With only one of the hadrons polarized, any nonzero spin asymmetry must involve at least one T-odd function. 
Restricting ourselves to hadrons with spin and 5, such functions do not show up in the transverse momentum 
integrated correlators $(x) and A(z). They do appear in the parametrization of the matrix elements involved in 
the decomposition of the transverse moments of the correlators. T-odd distribution functions only appear in the 
gluonic pole matrix element tt$Gi while T-odd fragmentation functions can appear in both the matrix elements Ag 
and 7rA(5 Using the parametrizations for these functions, one can calculate ei N -T,g(xi,x 2> Zi, z 2 ,y) and find 
the expression for the weighted cross section using Eq. I|22b|) . Considering only the (anti)quark contributions in 
p^+p — > ir+ir+X the resulting cross section is explicitly given in appendix [D] including in each term explicitly the 
factor Cq^ between braces { • }. 

The results can most conveniently be expressed as a folding of distribution and fragmentation functions, now 
including one T-odd function and a gluonic pole cross section 



( i sm(6(/>) da[hih 2 ] ) 

= dx 1± _ dx 2 ± drjx dr\ 2 cos(</>f ) / — - 
2?r J x± 



1/1 E rt^^fT^) I Df( Zl )Df{z 2 ) (30a) 



x±^ ^ 11 v ' 1 v ' 2 dt 
91929394 

M 2 , s,^_L(l)^ , S "^^9j39 2 ^9394 n 93^„ N n 94 , 



±= ]T hf(x 1 )h^ 1 \x 2 )^ ?1 2^ 3 ' 4 Df{z l )Df{z 2 ) 



XX , 

9192 93 9-1 



(30b) 
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E ^f(^)/f te^ ^V^ H^ {1 \zx)Df{z 2 ) +{K 1 ~K 2 ) (30c) 

v 9192 93 94 



9192 93 91 



(30d) 



where the summations run over all quark and antiquark flavors and the angle <fif is defined by <f>f = 0i— </> s . All 
non- vanishing partonic scattering cross sections and gluonic pole cross sections are functions of y or, equivalently, 
r\\—r\2 and those that contribute to hadronic pion production are listed in appendix El 

We note the occurrence of one T-odd function in each of the terms in Eq. i|3U|) . the functions f^ 1 \x) and hi\x) 
coming from the gluonic pole matrix element tt&g, the function (z) coming from the link-independent correlator 

Ag and the function H± (z) coming from the gluonic pole matrix element ttAq. We would, once more, like to 
emphasize that for fragmentation both Ag and ttAq contain T-odd functions that could contribute to the Collins 
effect. In Rcfs |^,|32 it is argued that the Collins effect is universal. This situation would occur if gluonic pole matrix 
elements in the case of fragmentation vanish, in which case the function vanishes and all T-odd effects come from 

the 'universal' function H^ x \ The latter function appears folded with ordinary partonic cross sections. However, in 
this paper we will allow for the gluonic pole matrix element for fragmentation and a nonvanishing function , 
one reason being that there is no full agreement between the different model calculations concerning the universality 
of the Collins effect |3J, y4( . 



B. Hadron-jet production in p'p scattering: p' +p —* 7r+Jet+X 



We only take into account (anti)quark scattering processes in the weighted scattering cross section for p^+p — > 
7r+Jet+X with the pion and the jet approximately back-to-back in the perpendicular plane. This cross section can be 
obtained from the more involved two-particle inclusive scattering cross section (IDlll by taking D\{zz) — 8{z2—X)5^ 2q = 
x± 5(x 2 ± — xj_)5 j2q and by letting all other fragmentation functions vanish. Here 5 J2Q is a delta function in flavor space, 
indicating that the jet j 2 is produced by quark q. The explicit expression using the diagrammatic approach is given 
in appendix [Dl This can be recast into a form involving distribution and fragmentation functions folded with gluonic 
pole cross sections 

{\$m(54>) do[hij 2 ]) 

= dx x ± dx 2 j_ drji dr\ 2 ^p- cos(</>f ) 

Z7T 



X 



E f«$ X \*i)f?{*2) I Df{z x ) (31a) 

x 2 ± V s * at 

91132 93 94 

+ E hf{ Xl )h^\x 2 ) S - ™>™ D?{ Z1 ) (31b) 

x 2 ± V s z at 

91929394 



^ £ hf{ Xl )fr{x 2 ) S - ?i^3±H^ {1) (z 1 ) (31c) 



X2 

91929394 



£ hf ( Xl )f? (a,)! dA<J ^^ H*t {1) M I (31d) 

X 2 ± y S Z clt J 

91929394 



C. Jet-jet production in p'p scattering: p' +p — > Jet+Jet+X 

We only take into account (anti)quark scattering processes in the weighted scattering cross section for p^+p — > 
Jet+Jet+A with approximately back-to-back jets in the perpendicular plane. As argued, in principle one can construct 
azimuthal asymmetries that give access to Eg (xi, X2,y), by weighting with the small momentum q". However, this 
requires accurate determination of the jet momenta k\ and k 2 . Here we only present the cross section obtained by 
weighting with sm(6<fi), which can be obtained from the more involved two-particle inclusive process (|D1(I by taking 
Df(zi) — S(zi~l)S : ' iq = x± S(xi±—x±)8 : > iq and by letting all other fragmentation functions vanish. Casting the result 
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from the diagrammatic approach (given explicitly in appendix [Dj) in the form of a folding, one obtains 
( | sin(^) da\jij 2 ] ) 

= dxu_ dx 2 j_ drji dr\ 2 cos(^)f) 5(xx±-x 2 ±) 
Ztt 

x{§ E f«±£H*i)f?i*2)i * Tqi y» q ' (32a) 

91929394 

+ S E ^(^)^i L(1) (^) ^ A ^ff^ 3g4 } (32b) 



q\qi 93 94 



V. SUMMARY AND CONCLUSIONS 



In this paper we have used the diagrammatic approach at tree-level to derive expressions for single transverse- 
spin asymmetries in 2-particle inclusive hadron-hadron collisions. The final states considered are hadron-hadron, 
hadron-jet and jet-jet, which are approximately back-to-back in the plane perpendicular to the incoming hadrons. 
The single spin asymmetries require the inclusion of transverse momentum dependence for the partons. We have 
assumed factorization to hold in this treatment of TMD effects although it is, at present, certainly not clear whether 
such a factorization holds for hadron-hadron scattering processes with explicitly TMD correlators. We have limited 
ourselves to the first transverse moments obtained by weighting linearly with the transverse momentum. These 
transverse moments show up in azimuthal asymmetries. 

While single-spin asymmetries generated by fragmentation processes, in which one can have T-odd fragmentation 
functions, are well-known, the single-spin asymmetries connected with initial state hadrons are more subtle. Within 
the diagrammatic approach T-odd effects for transverse momentum dependent distribution functions are attributed 
to the structure of the integration path in the gauge- link. This path depends on the specific hard process in which the 
correlator is used, explaining, for instance, the appearance of the Sivers function /{y 1 with opposite signs in SIDIS 
and DY |j, l35t l3q | . In the transverse moments of quark and antiquark correlators the effect of the gauge- link appears 
via the gluonic pole matrix element, which in the case of distributions is a T-odd matrix element giving rise to single 
spin asymmetries [53, • In this paper we show how the effects of the gauge-link appear as factors Cq , which 
determine the strengths with which the gluonic pole matrix elements occur. This is a generalization of the factors ±1 
appearing in SIDIS and DY. The fact that these strengths are determined by the hard parts makes it convenient to 
absorb them in so-called gluonic pole cross sections. Just as the transverse momentum averaged cross sections can, in 
leading order, be written as a folding of universal distribution and fragmentation functions and a hard partonic cross 
section, the single spin asymmetries can be written as a folding of universal distribution and fragmentation functions 
(involving one T-odd function) and a gluonic pole cross section. 

In our approach we allow for two possible mechanisms to produce single spin asymmetries in the case of fragmen- 
tation. This implies that in the two matrix elements in which the transverse moments can be decomposed, i.e. the 
link-independent part Ag and the gluonic pole matrix element tt Aq, one has both T-even and T-odd effects. For the 
Collins effect in fragmentation, it leads to two independent functions and H^^\ the latter appearing in the 

paramctrization of the gluonic pole matrix element. Having different linear combinations of these functions in SIDIS 
and electron-positron annihilation spoils the comparison of the Collins effect in these processes. In hadron-hadron 
collisions we find other linear combinations of the two functions. If fragmentation functions are universal, as is argued 
in Refs. 0,H2, the tilde function H^ 1 ' (and the gluonic pole matrix element for fragmentation) vanishes. In that 
case only the contribution from remains. 

Our results, including the strengths of the gluonic pole matrix elements differ from those of earlier calculations in 
which the effects of the gauge-links have been omitted. However, these effects can easily be incorporated by using the 
gluonic pole cross sections instead of the normal hard partonic cross sections. 

We have restricted ourselves to a particular single spin asymmetry in hadron-hadron scattering where the asymmetry 
arises from the deviation from the back-to-back appearance of the produced hadrons/jets in the perpendicular plane. 
This situation was discussed in Ref. ^} without inclusion of the effects of gauge-links. Although experimentally 
more challenging, the 2-particle inclusive case is easier to analyze than the 1-particle inclusive case, where large 
single spin asymmetries are observed, but where subleading transverse momentum averaged T-odd fragmentation 
functions |3?1 ] will also contribute. In principle the diagrammatic approach allows for inclusion of these contributions. 
Furthermore, the methods used in this paper to include the T-odd, transverse momentum dependent effects in 
(anti)quark contributions, which are crucial to treat single spin asymmetries in hadron-hadron scattering, can be 
extended to include the gluonic contributions as well as to treat various T-even double spin asymmetries. 
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FIG. 4: Gauge-links entering in the correlator for the lower-left incoming quark for a hard 
two-fermion scattering process without exchange of charge. Top: quark-quark scattering; 
bottom: quark-antiquark scattering. 



APPENDIX A: QUARK CORRELATORS AND GAUGE-LINKS 

The starting point for the structure of the hadron^quark transition is the quark correlator = P, S) |3Sll3£ 



d 4 t 



e l P<(P,S\^J0)UO\P,S) 



Similarly, one has for the quark— »hadron transition the fragmentation correlator A(fc) = A(fc; K, S), 



d 4 £ 



\MO<a h ^(0)\0) 



with 



d 3 P 



x 



4 a h ^ x (2,t)^2E a 



\Px\K, S h ){P x ;K,S h 



(Al) 



(A2) 



(A3) 



In the description of hard scattering processes we need the quark correlator and the fragmentation correlator integrated 
over, at least, the partonic momentum component p-P. This leaves the TMD correlator 



Hx, Pt ) = J d( P -p) . 

Integrating the TMD correlator over or weighing it with the transverse momentum p T , we obtain 
$(x) = J d 2 p T $(x,p T ) , 
$g(a:) = J d 2 p T p% $(x,p T ) . 



(A4) 

(A5a) 
(A5b) 



One finds similar expressions for the fragmentation correlator. Analogous to the above one can write down the 
antiquark correlator $ describing the hadron— >antiquark transition and the antiquark fragmentation correlator A 
describing the antiquark^hadron transition. 

To obtain properly gauge invariant correlators, gauge-links connecting the parton fields in the matrix elements are 
needed. The general structure of the gauge-links is U\ c ^ (0, £)=V exp[ — ig J c A(z)-dz], where the integration path 
C runs from to £. Here A is the gauge field and V is the path-ordering operator. The integration paths can be 
calculated by rcsumming all collinear gluon interactions between the soft and hard parts. Consequently, for the TMD 
correlators they depend on the process in which they occur. The gauge-links appearing in the quark correlator in 
a two-fermion hard scattering process with uncharged boson exchange, such as in QED, are readily calculated by 
considering the flow of the fermion lines [2(| . The results from that reference are given in Fig. 0] Explicitly, we 
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encounter the link structures 



^[(0-,0 T ),(±cx)-,O r )]^[T±oo-:0 T ) : (±oo-,cx> T )] C/ ['[±oo-,oo T ) : (±oo-,$ T )] t/ [(±oo- 



,£t)>(«-,€t)] 



which are build up from the gauge-links along straight lines 

fb 



HI 



dz n ■ A(z) 



and u [a,b] = P ex P 



dz T ■ A T (z) 



(A6) 
(A7) 

(A8) 



The gauge-links in the scattering of two colored fermions in QCD can be obtained from those in Fig. 01 by accounting 
for the flow of color charge, using well-known QCD rules for color flow such as 



= Tf 



1 



For example, the tt*-channel of quark-quark scattering can be decomposed in this way giving 



Ti 




1 

~N~r. 



1 

N 2 



(A9) 



(A10) 



The gauge-link of this diagram can be obtained by replacing each diagram on the r.h.s. with the corresponding QED 
gauge- link as given by Fig. 21 and factoring out the overall color factor of the QCD diagram. The overall color factor 
of the gluon exchange diagram on the l.h.s. is (Tr[t a t b ]) 2 —Tp(N 2 ~ 1), which can also be obtained by tracing the color 
flow in all diagrams on the r.h.s. This color factor does not enter in the gauge-link, but in the evaluation of the 
diagram itself and is included in the hard amplitudes that will be used in the calculations in appendix [Dj Accounting 
for the additional factors Tr(l) = N c that are obtained for each color loop, one obtains the gauge-link 



T 2 (N 2 



l)xU%p=Tl[Nl 



■qq -^ L -c j^^ 

The other diagrams can be calculated analogously. For quark-quark scattering we obtain: 



tsiqq qq 



N 2 
N, 



L_{ (iVc 2 + 1) - } 



;/[*"*] _ Jj[ut*] _ 
99 99 jy2 



2N r 



Nr. 



Nr 



US 



and for quark-antiquark scattering: 



U 



N 2 - 



j/tt*} _ 



-A 



Tr(U 



N c 
Pltl 



(Nc 



U. 



St*] 



. ?A ts ' ] - Nc 



i{ 



2)WH}, 

Nr. J 



(Alia) 
(Allb) 

(Al2a) 
(A12b) 



99 -99 N^-ll"" N t 

These are the gauge-link operators that enter between the quark fields in the correlator of the incoming quark: 



(Al2c) 



2tt (2tt) : 



e*-e(P ) s\i>(PMo,t)ii>@\p,S) 



(A13) 



The gauge-links that enter in the quark-fragmentation correlators are the time-reversed ones as compared to those in 
the quark-correlators. That is, aW' + ' in the quark-correlator corresponds to aM' in the fragmentation correlator and 
a to a etc. The gauge-links that enter in the antiquark-correlators $ and A are the hermitian conjugates 

of the gauge-links in the quark-correlators $ and A of the corresponding diagrams. 

We note that for the fragmentation correlators the gauge-links are all split up in parts, parts connecting to the field 
at £ and others to the field at 0. Taking as an example the quark fragmentation correlators with the gauge-links W'^ , 
one has (compare with Eq. l|A6jl ) 

d(i-K) d 2 ^ 



A^(z,k T ;K, S) = J 



2tt (2tt) 



U[(±00 - , 0O T ) , (±00 - 4 T )] U [( ± C 



x ala h <ip{0)U [(0 - OT 



,«t),(€-.€t)1 

V7 



(A14) 



),(±oo-,0 T )] Ly [(±oo-,0 T ),(±oo-,oo T )] 
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APPENDIX B: CONSEQUENCES OF GAUGE-LINKS FOR DISTRIBUTION FUNCTIONS 

The gauge-link has important consequences for the parametrizations of the correlator in T-even and T-odd functions. 
We start with the link structures enumerated in Fig.Q] The TMD correlators are link dependent. We write for the 
correlators with gauge-links W^l, $P+] for U^U^ + \ $[( D )+1 for ^-Tt(U^)U^ and $K nt )+] for Tr(U™)U^+l 
We will also need the transverse momentum integrated correlators $2, (a;) an d &q(x, x—x') 

*£(*) = / e fa ^ p )(P,5|^(Q)f/- S] ^(0V(e)|P,5')J LC , (Bl) 

$< G {x,x-x>) = J d{ ^p d{V 2 p e^-^ p V^ p )(P,^^ , (B2) 

which are set on the lightcone (LC) where £ ■ n = £ T = and i] ■ n — r\ T = 0. We have also used the shorthand notation 
G na = g fl „G tia n 1 ' for the held strength tensor. In terms of these the weighted correlators can be written as 

* [ B ]a (x)=*%{x)- I dx' - 7 ^$%(x,x-x') = $%(x)±w$%(x,x) , (B3a) 

$ l g +]a (x) = $£(ss) - J dx' { - ^(x') } $g( x > x ~ x ') = ®d( x ) + ^g( x , x ) , ( B 3b) 

*a D)+1 "(^) = *L (Dt)+lQ ( a; ) = $ SO) - / dx> -y^- a?-*') = ^(i) + 7r*£(z, a?) , (B3c) 

J X 16 

where $3 without link index refers to 

$g(a;) = - / dx' P-L ^fox-a/) . (B4) 



The decomposition in Eq. 1|B3(I is useful because time reversal symmetry implies that the correlator $g only contains 
T-even functions, while <3?g only contains T-odd functions. 

The correlators encountered in p^p — > irnX are readily obtained from the results above and can also be decomposed 
in terms of <£>a and For instance, for the £t*-channel in qq scattering we get from Eq. (|AllaJ) 



= (^c + 1) - 2 } + (N? + 1) - 6 } TrSgfoa;) . 

The other correlators can be calculated analogously. For qq scattering we obtain: 

= = + ^l^g^x) , (B5a) 

= * [ f ]a (x) = - ^f±! 7rd>£(z, s) , (B5b) 

and from Eq. I|A12|) we get for qq scattering 

4 SS * ] » = * [ f ]a ( x ) = * [ ?* ]a ( x ) = ^( x ) + ^f^^ G ( x , x ) > (B6a) 

4 tn » - - ^f^TrSgfoa:) . (B6b) 
The integrated quark correlator <j>(cc) is parametrized as follows in terms of quark distribution functions pil I40L 

mm 

$u(a;;P) = \ h{x)f , (B7a) 

$ L (o:;P) = i^ 5l (x)75f , (B7b) 

$ T {x;P) = \h l {x)~ f5 \[$ T ,f], (B7c) 
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Uqq 


7 /[*«*] 
Uqq 


7 /["'*] 
Uqq 


Uqq 


7 /I ss *] 
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7 /[*"*] 
"99 




$[«*] 
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jjj[tU*] 
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<j>[ss*] 
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2v|-l 


ivf-1 
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TABLE I: The basic gauge-links (upper table) and the gauge-links in specific hard scattering (qq and gg) 
diagrams (lower table), the notations used for the correlators and the strengths Cg of the gluonic pole 
contribution tt<&g- 



where 

e ^ = -^-e Pn ^ : and S = S L ±- P - S L7 ^§- n + S T , (B8) 
P-n M 2 P-n 

with Sl+S^, = — 1. The indices U, L and T refer to unpolarized, longitudinally and transversely polarized hadrons, 
respectively. For the T-even transverse momentum weighted correlator $9 (a;) and the T-odd gluonic pole tt^>g{x,x) 
one has the parametrizations 

= , (TrSgy^P) = |M^«(x)i[f )7 a ] , (B9a) 

(*g) £ (x;P) = ^M^^hgitf^l , {ir^ G ) L {x-P) = , (B9b) 
(*g) r (z;P) = iMS« 5 W(x) 75 f , (7r$g) T ^;P) = iMef/^^f . (B9c) 

From the parametrizations given above and using the decomposition in Eq. (|B5a|) . we find that the T-odd distri- 
bution functions /{y 1 ' and hf^ appear with a multiplicative prefactor Cq 1 —(N^—5)/(N^—l) in the contribution 

\u [±] ] 

corresponding to the tt*-channel in gq-scattering. This is the appropriate generalization of the factors C G = ±1 
occurring in SIDIS and Drell-Yan scattering (as explained in section IIIII . Similarly, the prefactors of the T-odd 
distribution functions appearing in the other scattering channels can be read of from Eq. (|B5|) for qq scattering and 
from Eq. i|B6() for qq scattering. These prefactors are summarized in TableQ] From the Eqs. (|B5|) and l|B6|) we also see 
that all the T-even distribution functions occur in hadron-hadron scattering in the same way as they do in SIDIS, i.e. 
with a prefactor +1. For antiquark distribution functions, which can be related to quark distributions in the negative 
x region, the same results as above apply. The antiquark distribution functions will be distinguished from their quark 
counterparts by an overline, e.g. fi(x), etc. 



APPENDIX C: CONSEQUENCES OF GAUGE-LINKS FOR FRAGMENTATION FUNCTIONS 

The discussion on the consequences of the gauge-links for fragmentation functions is a little bit more involved than 
for distribution functions, due to the presence of the hadronic states \K, X) in the definition of the correlators. These 
are oui-states, preventing the use of time-reversal symmetry to constrain the parametrization. 

All collinear interactions between the soft and hard parts result in the quark-fragmentation correlator £^~\k) in 
SIDIS and the correlator AM(fc) in electron-positron annihilation (see equation 1A14(1 ). The transverse-momentum 
integrated fragmentation correlators in these two processes are 

A [ ±](z) = J d(^l e ^ K y^0\U [M] ma{a h mU^ ±oc] |0)J LC . (CI) 

Although not immediately evident, it is not hard to see that, since there are only gauge-links along the n/j-direction, 
the two correlators are identical: AW(z) = AH(z) = A(z). 

In analogy to the previous appendix we define a correlator A^, and a gluonic-pole matrix element A G 

A£(z) = J ^^e i (W'(0|^ i gii? a (OV'(04a^(0)^ [ o I c]|0)J LO) (C2) 
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A Q fi l—L) = f d ^' RS > d (T K ) i(t-K)/z A\{yK)-(i-K)}/z' 

GV *' Z z " J 2tt 2tt 

x (0| gG n - a ( V ) m a{a h ?(0) |0) J LC , (C3) 

with G n,lQ = g^G^n^ (and £ an arbitrary point). It can be shown that in terms of these the weighted correlators 
can be written as 

&£ ]a (z) = A a D (z) - [ d(jr) ^— Ag(±, i-i) = Ag(*) ± 7rA£(±, i) , (C4a) 

J 2 'T JC 

A [ d Di]a (z) = A" (z) - | d(±) { ^— + 2^) | Ag(I, i-i) = Ag(z) - 37 rAg(I I) , (C4b) 
Al- (D)ltt (z) = Al- (Dt)lQ (z) = AS (z) - [ d(jr) *- Ag(i, i-i) = Ag(z) - 7rAg(i, I) , (C4c) 
where A3 without link index refers to 

Ag(*) = A£(*) - / dz'- 1 Ag(i, i-i) . (C5) 



As stated at the end of appendix E] the gauge- links in the fragmentation correlators in p^p — > 7T7tX are obtained 
from IjAllfl and (|A12I) by time-reversal. We, then, find the following quark-fragmentation correlator for the tt*-channel 
in quark-quark scattering (cf. IjAllal) ): 

Af ] » = j^-^ (N? + l)A^(z) -2A^\z) } 

= ^ L T {(^ + l)-2}A^)-^ T {(JV c 2 + l)-6}^(I,i). 



The other quark-fragmentation correlators can be calculated analogously. For quark-quark scattering we obtain: 

A f l^^M'W^^-l^^^i) , (C6a) 
A^>(z) = Af >(*) = Ag(z) + 7rAS(i, I) , (C6b) 

and for quark-antiquark scattering 



Ajf >(z) = Af lQ (z) = A^ lQ (z) = Ag(z) - ^±1^(1, i) , (C7a. 



Af l«(z) = Ag(z) + ^-|7rAg(i,i) . (C7b) 

The integrated fragmentation correlator A(z) is parametrized as follows 0] 

zA u {z;K)=D l {z)f, (C8a) 

zA L (z;if) = S t Gi(z) 75 # , (C8b) 

zA T (z;if) = H 1 {z)^\[$ T ,$] , (C8c) 



with 



' Kn hl xu A o - o _L Ef _ Q. Mfe 



c 



and 5 = 5,,— X- S L —?—n h + S T . (C9) 
K -n h M h 2 ii -n/j 

The functions in these parametrizations are called quark fragmentation functions. Due to the internal soft interactions 
in the final-state hadron the correlators Aq and ttAq both contain T-even and T-odd parts @. Correspondingly, 
they have very similar parametrizations in terms of fragmentation functions. We will distinguish the fragmentation 
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7 /["**] 
"99 


7 /[ts*] 
"99 




A [«*] 


A [««*] 


A [tu*] 


A [«t*] 


A [tt*] 


A ( ss *] 


A [«t*] 


A [ts*] 


°G 


1 
2 


1 
2 


3 
2 


3 
2 


3 

4 


5 

~~ 4 


5 
4 


5 
4 



TABLE II: The basic gauge-links (upper table) and the gauge-links in specific hard scattering (qq and 
qq) diagrams (lower table), the notations used for the correlators and the strengths Cg of the gluonic pole 
contribution nAc with N c =3. 

functions in these two correlators by adding a tilde to the fragmentation functions appearing in the parametrization 
of the gluonic pole. That is, parametrizing the correlator Ag as follows 



z{A%) u {z-K) = M h iHl (1 \z)\[$ n a ] , (ClOa) 

z {A%) L {z;K) = S ll M h H$\z)7t\[$,'f l ] , (ClOb) 

z (A$) T (z;K) = M h [s«G ( £(z) lb $ - ef T D^ l \z)%} , (ClOc) 

the parametrization of the gluonic pole is written as 

z(^Ag) a (I,±;X) = M h iH± {1 \z)\[$ n a ] , (Clla) 

z ^A a G ) L {\,\;K) = S L M h H^ 1} (z) l5 ^ 7 a ] , (Cllb) 

z^A%) T {\,\ ] K) = M h [s^G^{z) 1 a-^f T D^ 1) {z)$) . (Cllc) 



The fragmentation functions appearing in these parametrizations contribute to azimuthal asymmetries in special 
combinations. For instance, using the decomposition in Eq. (|C6a|l we find that the Collins effect contributed by 

the tt*-channel for qq scattering is ifp 1 ' — jftzr H^ 1 ' . Similarly, the other partonic channels contribute particular 
combinations of fragmentation functions. Which combination of fragmentation functions one should take for a certain 
process can be read of directly from the decompositions in Eq. lC6p and l|C7|l . That is, if we let FF(z) denote a generic 
fragmentation function appearing in the parametrizations in Eq. I|C10|) and IjCllll . then this fragmentation function 

will appear in the expressions for azimuthal asymmetries in the combination FF(z)— Cq^ FF(z). In particular, we 

see that the tilde fragmentation functions always appear with the (process dependent) prefactors cff 1 summarized in 
Table UTI while the fragmentation functions without a tilde always occur with a simple prefactor +1. If the gluonic pole 
matrix elements ttAq vanish, then so do all the tilde functions. In that case fragmentation is completely described 
by the universal functions appearing in the parametrization of Ag. Notably, the Collins effect is always given by the 
term H^ (1 \z). 

For antiquark- fragmentation functions, which can be related to the quark-fragmentation functions in the negative 
z region, the same results as above apply. 

APPENDIX D: RESULTS IN THE DIAGRAMMATIC APPROACH 

In the expressions given below y is given by Eq. J7J and s is s = x\ s cosh 2 [5(771— %)] — x\s/Ay{l—y). The 
summations run over all quark flavors, including the case that q' = q (where applicable). Similarly, the S Jiq are delta 
functions in flavor space, indicating that the jet ji is produced by quark q. We have written the factors Cq^ for 
the T-odd distribution functions between braces {■}. For the Collins functions we have written the combinations 
jjM 1 ) —Cq H^ 1 ^ 1 between braces. The factors Cq are taken from Table |U for the distribution functions and from 
Table ITT1 for the fragmentation functions. 



+p — > TT+ir+X 

Considering only the (anti)quark contributions in +p — > 7T+7T+X the resulting cross section is given by 

( 5 sin(5<^) dcr[/ii/i 2 ] ) 

= rfa;i_L dx 2 j_ ^p- cos(0i-0 s ) [ 

27T J x± 



x 



iZ^i f „„ (i-y) 2 + i 

9 s 2 



{ M i (1 y + 1 EU}/^r 1) (^i)/?'(^)£>?(zi)£>? , (^) 
" 9,9' 

+Mx ((1 - yf + y 2 ) £{f} f^ 1 \x 1 )fl(x 2 )Di(z 1 )Dl(z 2 ) 
9,9' 

+M 1 (1 ~ V X + 1 ]T{-f l/^Cxr)/!' (x 2 )^(. 1 )^'(z 2 ) 



y 9,9' 
\2 



+M 1 |iL^L^ { | }/ ^a) (xi)/ - Kx2) ^ (zi)5?(z2) 

d y q 

-M 2 2y(l y) £{§} h\{ Xl )h^ (1) (x 2 )d( ( Zl )D\(z 2 ) 



g,q' 



-M 2 -(1-y) ]T{f } K{ Xl )h^ \x 2 )Dl{ Zl )D\{z 2 ) 



~M 2 ^ ^^{— §} h\{xi)h q ^ 1 \x 2 )D\{zi)D\{z 2 ) 

9 



-M, 1 2^^/ i ?(x 1 )/Y( a;2 ){^ (1) (^)+!^ (1) (^)}5f(^) 

^ 9,9' 

-m^I^E^^^^i^^i)-!^ 1 '^)}^^) 

J ^ 9 

+mJ^/ 1 «( Ii )/?(, 2 ){^ (1 »(. 1 )45^( 1 »(, 1 )}c?(, 2 ) 



3 

9 



+ ( quarks <-> antiquarks ) | + ( K\ <-> i^ 2 ) 
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p + p -> 7T + Jet + X 

Considering only the (anti)quark contributions in +p — > 7r+Jet+X the resulting cross section is given by 
( 5 sin(^) da[hi j 2 ] ) 



dxij_ dx 2 j_ ^771 d?7 2 ~r ~~ cos(0i— <j> s ) 



s 



2ir 2x 2 ±^/s 
„2 r n _ .^2 



x 



%^ { M 1 (1 ~ + 1 ^{1} rt£\xi)f< {x*)]*M (D2a) 

9,9' 

+ M i7f^EU}/ 9 ^ 1) (^)/f'(^)^ 1 '(^i)^« (D2b) 
1 yj 9,9' 

-M^-^^^-l}/^ 1 )^)/! 9 ^)^?^)^ (D2c) 

+Mi ((1 - y) 2 + y 2 ) ]T{f } rtP(xi)n(x 2 )Df ( Zl ) 5™' (D2d) 

+Mi ((1 - y) 2 +y 2 ) ^{l}/^^!)/?^)^'^)^^' (D2e) 

9,9' 

+M! (1 'g 2 + 1 EM > f q ^\^)H{x2)D\{z{) 5™' (D2f) 



^ 9,9' 



+M 1 ^^^{-l}r^ 1) (x 1 )n(x 2 )Di(z 1 )5^ (D2g) 



(i-y) 



9,9' 
\2 



+Ma liL^-^l}/^ 1 )^)/!^)^^)^ 9 " (D2h) 

^ V 9 

+M 1 ^ T ^^{|}/^ 1) (x 1 )/l(x 2 )5?(z 1 )^ 9 (D2i) 

^ 9 

-M 2 22/(1 - y) £{§} Z^)^ 1 ^)^' (*) (D2j) 

9,9' 

-M 2 2y(l - y) ^{f } Z^)^ 1 ^)^' (^) (D2k) 



-M 2 |(l-y) E{|} hl( Xl )h"t {1) (x 2 )Dl( Zl ) 5™ (D21) 

9 

-^^{l}^^)^^)^^)^ 9 (D2m) 

9 

-M 2 | hl( Xl )h*tW (x 2 )D\{ Zl ) 5™ (D2n) 

9 

-M.^^^^^O/f'^)!^ 1 )^)-^^^)}^' (D2o) 

9,9' 

+^i|^E^(^)/i'(^){^ (1) (^i)+i^ (1) (^)}^ ( D2 p) 

^ 9 

-M hl 2^^hl(x 1 )fi(x 2 ){H^ 1 \z 1 )+^H^ 1 \z 1 )}S^' (D2q) 

9,9' 



19 



-M^l^^^^)/?^)!^ 1 '^)-!^ 1 '^)}^' (D2r) 



3 y 



+M hl lY / K( x i)n(^) {H^ {1 \z 1 )-^f {1 \z 1 )} (D2s) 



3 

9 



- ( quarks <-> antiquarks ) 



p + p -> Jet + Jet + X 

Considering only the (anti)quark contributions in +p — ► Jet+Jet+X the resulting cross section is given by 
(isin(^)dCT[jij 2 ]) 

= dx 1± dx 2 ± d-qi drj 2 ^p- cos(0i-0 s ) S(x 1 _ L -x 2 ±) 



x 



47r «s f (l-2/) 2 + l 



9 s 2 



| Ml (1 ^ 2 + 1 ^{1} f*#\xi)tf (x 2 ) 5^8™' (D3a) 

" Ml 5^(1^) £H> (D3b) 
+M X ((1 - y) 2 +y 2 ) £{f}/ 9 # 1) (*i)/fc 2 )<P 9 >«' (D3c) 

9,9' 

^ 9,9' 

+M X £{f } f^\ Xl )n(x2) (D3e) 

^ ^ 9 

-M a 2y(l - y) ]T{| } ^i)/^ 1 ^) ^V 2 *' ( D 3f) 



9>9' 

"M 2 |(1 - y) £{f } MCasi)^^ 1 '^) 8^8^ (D3g) 



-M a i £{-§} fc?0»i)fc*i W (*a) 5 il9 <F 29 (D3h) 

9 

+ ( quarks antiquarks ) | + ( Jet x <-> Jet 2 ) 

APPENDIX E: PARTONIC CROSS SECTIONS 

In this appendix we enumerate all the (anti)quark scattering cross sections (taken from |44|) that are needed in this 
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Quark-quark scattering 

The unpolarized quark-quark scattering cross sections are given by 



dcfqq'-^-qq' 


Aira 2 


s 2 


+ 


u 2 


di 


9 s 2 




i 2 


3 


do ' qqi —+q'q 


Aira 2 


s 2 


+ 


P 


di 


9 s 2 




U 2 




do qq-^-qq 


_ dcTqq'-> 


qq 1 


+ 


d(7qqi —tq'q 


di 


di 






di 



d&qq—tqq 

di 



where da 1 represents the interference terms 

di ~ 27 s 2 tii ' 
The polarized quark-quark scattering cross sections are 



dAd-qi q i_ qq _ 8 



dt 

dA<T q T g> g> 

di 



27 s 2 ' 



■snaZ us 



9 s 2 P ' 

dA6- oTo ^ T dAa q i q ,^ q i q , dAa 1 



qlq^q'q 



dt dt 
with the interference term 



dt 



di 27 s 2 t ' 

The modified cross sections are 



d&gqq'^qqi _ N 2 — 5 da qq i ^ qq i 



dt 



N?-l di 



d(Jgqq'—>qiq N 2 5 d&qq' -^qi q 



dt 



n 2 -i di 



dOgqq^qq _ N 2 ~5 



dt 

d^q^gq^qq 



N?-l 



da„ 



do 'qq' -^q' , 



dt 



dt 



2 iV c 2 +3 da\q^ qq 

n?—i di 



N 2 +3 dAd q ,q^ qq 



di N 2 -l di 

dAa q rg,^gqiq, _ N 2 -5 dAa q T q ,^ q t q , 
di 

dAa q r g^ggt q 



n 2 -i di 

N 2 -5 dAa q , q ,^ qU , iV 2 +3 dAa l qU ^ q 



di N 2 -l di N 2 -l di 

The partonic cross sections can be regarded as functions of the variable y defined in Q) through 

i 



= -y 



= -(i-y) 



9 

X 1 S 



4y(i - y) 



(Ela) 
(Elb) 
(Elc) 

(E2) 

(E3a) 
(E3b) 
(E3c) 

(E4) 

(E5a) 
(E5b) 
(E5c) 
(E5d) 
(E5e) 
(E5f) 

(E6) 



All other non-vanishing quark-quark and antiquark-antiquark scattering cross sections that contribute to hadronic 
pion production can be obtained from the above from symmetry considerations. 
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Quark-antiquark scattering 



The unpolarized quark-antiquark scattering cross sections are given by 

d<T q q'^ q q' 4:TTa 2 S 2 + U 2 



dt 



d& a 



9 s 2 t 2 
4wal i 2 + i 



dt 



9 s 2 



d& a 



d& a 



dt dt 
with the interference term 



+ ■ 



dc Tgq^q'q' ^ ^ qq^qq 



dt 



dt 



da. 



qq^qq 



Ana 2 u 2 



di 27 s 2 ts 

The polarized quark-antiquark scattering cross sections are 



dA<7 gT? T 


-^q'q 


8wa 2 s 


iu 




di 




9 s 2 






dAtJqt q1 


^qq 


dAa q tqt 


->g'<f 


dAa 1 ,-, - 

qTqT^qq 


di 




di 




di 


dAa q Tq,_ 




8wa 2 s 


us 




di 




9 s 2 






dAa q rq_ 


>q1q 


dAaq\-q,^ q rq, 


dAa 1 ,- ,- 
q<q->q< q 


di 




di 




di 


dAa q i q _ 


>qtq 


8ira 2 s 






di 




27s 2 







with the interference terms 



dAa 1 ,-, - 
di 



8ira 2 u 
' 27 s 2 1 



dAa 1 ,- 

q~q^q~q 

di 



8na 2 u 
27s 2 1 ■ 



The modified cross sections are 



da- n 



N 2 -3 da a 



dt 

d&qqq^q'q' 



N 2 -l dt 
N 2 + l da q q^q,q, 



dt 



da 



gqq^qq 

di 

dAa q , w _ 



N 2 -l di 
N 2 —3 da, 



qq'^qq' 



N 2 -l di 

*q'q' _ N c+ l dAaqTq-r 



N 2 +l 
N 2 -l 

->q'9' 



da. 



qq^q q 
di 



da, 



qq^qq 
di 



dt 

dAa q , w 



N?-l 



dt 



+qq 



N 2 + l dAa oT(5 T_» 



TO 



dt 

dAa q ]q,^gq1q, 

di 

dA(T q 1 q^gql q 

di 

dAa qrq ^~ qrq 
di 



N 2 -l 



dt 



N 2 -3 dAa q rq,^qiq, 



N 2 -l dt 
N?-3 dAa. 



q'q'^q'q' 



N 2 -l dt 

N 2 +l dAagyg. 
N 2 -l di 



+ 



/V 2 +1 dAa 1 ,- ,- 

-iv c -p± q \ q ^ q \ q 



N?-l 



dt 



All other non-vanishing quark-antiquark scattering cross sections that contribute to hadronic pion 
obtained from the above from symmetry considerations. 



(E7a) 
(E7b) 
(E7c) 

(E8) 

(E9a) 
(E9b) 
(E9c) 
(E9d) 
(E9e) 

(E10) 

(Ella) 
(Ellb) 
(Ellc) 
(Elld) 
(Elle) 
(Ellf) 

(Ellg) 
(Ellh) 
production can be 
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